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Abstract 

Thermal emission rates and mean free paths of photons in a color-flavor locked (CFL) 
phase of quark matter at high densities and moderate temperatures are evaluated. Our 
calculations are based on a low-energy effective theory for CFL matter describing Gold- 
stone boson excitations and their electromagnetic as well as strong interactions. In- 
medium coupling strengths of vector mesons are estimated to be smaller than in vacuum. 
As a consequence of in-medium modified pion dispersion relations, novel processes such 
as Tr~^TT~ — > 7 and 7 vr^vr" become possible. The total photon emissivity is found to be 
very large, exceeding contributions from thermal e^e~ annihilation above temperatures 
of about 5 MeV. At the same time, the corresponding mean free paths become very 
small. Our results imply that the photon flux from the surface of a (hypothetical) CFL 
star in its early hot stages saturates the black-body limit. Schematic estimates for the 
early thermal evolution of the star are also presented. 



1 Introduction 

The discovery of asymptotic freedom, leading to the formulation of quantum chromodynamics 
(QCD) as the theory of strong interactions, was soon followed by the suggestion that matter at 
sufficiently high densities consists of a deconfined phase of quarks [T]. Only shortly afterwards 
it was pointed out [2] that the true ground state of cold dense quark matter exhibits color 
superconductivity (CSC), characterized by diquark condensation with an estimated energy 
gap A of the order of 1 MeV between the highest occupied and the unperturbed quark state 
at the Fermi surface. Since this magnitude of the gap is rather small for phenomenological 
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applications, CSC subsequently received little attention. The situation changed when reinves- 
tigations OHj using nonperturbative forces {e.g., instanton- induced) showed that the gap can 
be substantially larger, A ~ 100 MeV for moderate quark chemical potentials, ~ 350 MeV. 
Similarly large values are obtained from estimates based on perturbative calculations at asymp- 
totically high densities E] • Thus, from the practical point of view, the existence of color 
superconductivity in compact stars has (re-) emerged as an exciting possibility. 

The detailed properties of CSC matter relevant to astrophysical applications depend on the 
interplay of the quark chemical potential, the q-q interaction strength, and the bare masses of 
the (light) quarks m, d and s. In particular, for Hq below the (constituent) strange quark mass, 
only u and d quarks are subject to BCS pairing. The corresponding phase is known as 2-flavor 
CSC (2SC). In the idealized case where the quark chemical potential is much larger than the 
strange quark mass (m^), the latter becomes negligible and all three flavors exhibit likewise 
pairing. The preferred symmetry (breaking) pattern in this phase corresponds to the so-called 
color-flavor locking (CFL) j7j, since the underlying diquark condensate is invariant only under 
simultaneous color and flavor transformations. In the present work, we will focus on the CFL 
phase (for a recent review and a more exhaustive list of references, cf. 0). 

The investigation of compact stellar objects largely relies on electroweak emission spectra 
over a broad range of frequencies. Thus, emissivities and mean free paths [mfp^s) of neutrinos 
and photons are important quantities to assess bulk properties of matter within compact stars. 
It is therefore vital to address the question in how far the existence of CSC inside compact 
stars affects their properties, such as neutrino and photon emission as well as thermal evolution. 
Pertinent studies have recently been carried out for CSC phases in various scenarios jm ilO| ITT]. 
Ref. [H] mostly addressed neutrino emission for CFL matter at temperatures in the sub-MeV 
range (~0.1 MeV or 10^ Kelvin) relevant to the long-term evolution of a compact star. In 
Ref. PU] the CFL neutrino problem has been investigated for a temperature range of up to 
30 MeV, corresponding to the hottest phases expected to occur during a supernova explosion 
or in the early stages of a proto-neutron star (see, e.g., Ref. ^21)- Photon emission has been 
calculated in Ref. [TT] in a weak-coupling expansion for (gapped) quarks. The objective of 
the present article is to evaluate both photon emissivities and m/p's in the strong coupling 
regime of CFL matter in the (tens of) MeV temperature region making use of an underlying 
effective theory description. The main processes of interest therefore involve scattering and 
decay processes of Goldstone bosons, taking into account their in-medium modified dispersion 
relations. For comparison, the effects of thermal electrons and positrons are assessed, which 
prevail towards lower temperatures. Furthermore, we incorporate generalized vector mesons, 
relevant at the high-temperature end of the CFL phase. As a lower limit, and in order to 
minimize the number of poorly known parameters, such as meson masses, we restrict the 
effective theory description to the SU{2) flavor sector, i.e., neglect contributions from mesons 
involving strange quarks (cf. Ref. JH] for a recent analysis in hot hadronic matter). We will 
comment on the effects of kaons and other strangeness-bearing mesons in the summary. 

The outline of the article is as follows. In Sect. El we will recall the basic elements of the 
effective low-energy Lagrangian which describes the electromagnetic and strong interactions 
of pseudoscalar Goldstone bosons. In Sect. IHl the evaluation of photon emission rates and 
mean free paths will be presented, including novel aspects specific to the CFL setting. Sect. El 
is devoted to an estimate of the photon-driven cooling behaviour in the early stages of the 
thermal evolution of a compact star which accommodates CFL matter throughout. Finally, 
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Sect. El contains our summary and conclusions. Explicit expressions for medium-modified 
Feynman rules and process amplitudes have been relegated to the Appendix. 

2 Effective Low Energy Lagrangian 

The effective low-energy theory of CFL diquark matter is analogous to the chiral theory in the 
QCD vacuum and has been introduced in Refs. ^3 El . Most notably, color- flavor locking 
induces the breaking of chiral and baryon number symmetries. Thus, an octet and singlet of 
pseudoscalar Goldstone bosons emerges which constitute the low-energy excitations and hence 
determine the features of CFL matter for energies below the gap, including electromagnetic and 
thermal properties. Although the (pseudo-) Goldstone bosons in CFL matter carry the same 
quantum numbers as in vacuum, differences in the chiral symmetry breaking pattern render 
their masses - arising form the explicit breaking due to finite current masses - qualitatively 
and quantitatively different. While in the vacuum case the parameters of the chiral Lagrangian 
have to be extracted from experiment, in the case of CFL matter a matching of the low-energy 
sector with the so-called high-density effective theory (HDET) allows for a determination of 
parameters in lowest order of the chiral and perturbative expansion, respectively. A detailed 
analysis for excitations below and above the gap revealed, e.g., that the mass hierarchy in the 
strange and nonstrange sector is (partially) inverted [ini HZl HE] ■ 

The main objective of the present work is an assessment of electromagnetic and strong 
processes involving Goldstone-boson and vector excitations. As a convenient approach towards 
this goal, we here adopt the hidden local symmetry (HLS) framework that, in the context of 
CFL matter, was suggested in Ref. and further elaborated in Ref. ^T]. It amounts to a 
treatment of vector mesons as dynamical, composite gauge bosons corresponding to a hidden 
local SU (3)c+L+_R symmetry.^ Originally, the HLS approach was introduced in the vacuum case 
in Ref. [20] (for a review see 21j). Here, we will briefly recall the main elements important for 
our work. 

We start with the effective chiral Lagrangian of (generalized) pseudoscalar Goldstone bosons 
in the CFL phase, which is equivalent to the non-linear sigma model known from the vacuum 
and reads [H 113 113 

C^^ = Y Tr [doUdoU^ - vl diU diU^] - c [det(M) Tr (M-^U) + h.c] , (1) 

with M = diag(m„, m^, m^) being the current quark mass matrix. The chiral field U is related 
to the Goldstone boson octet tt = 7r"A" by 

U = e'^l^^ , (2) 

where the Gell-Mann flavor matrices are normalized as Tr[A", A^] = 2 8°^ and is the (in- 
medium) pion decay constant. The latter, as well as the constant c in Eq. (Q), can be obtained 
from the matching procedure, referred to above, to leading order in «<, [TE] . 

2 ^ 21-8 1n2 ^\ ^ 3A2 

1 8 27r2 ' 27r2 ■ 

^Strictly speaking, hidden local symmetry only holds in the zero size approximation, neglecting the finite 
size of condensed diquarks, cf. Ref. ^H] for a discussion of this issue. 
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The explicit expressions for the pion masses then take the form jTH] 



ml -ml 
= - + 



3A2 , , 11/2 

rriu + md) m^ 



(4) 



As Lorentz invariance is broken due to the rest frame of the medium, the pion decay constant 
(defined via the weak pion decay matrix element) generally splits into temporal and spatial 
components j22|, denoted by fx and /s, respectively. They are related to the pion velocity 

(found to be = l/-\/3, coinciding with the speed of sound in relativistic fluids 
through fs = vlfr, with /t = f.„. Consequently, the pion velocity also appears in the effective 
Lagrangian ((H), and induces a modified dispersion relation for the Goldstone bosons, which 
follows from the corresponding kinetic term as 



E = ^ml^ + p2 . (5) 

The hidden local symmetry of the low-energy effective theory, Eq. (^, can be made explicit 
by decomposing the left and right pion fields 'nL,R according to 

UL,R = e^^^^-'f- = eciL,R. (6) 

where and ^l,r are SU{?>) color and chiral gauge fields, respectively. In the following, we 
will fix the local SU{?>)c+l+r gauge by = ^L,i? = ^5 which identifies the left and right pion 
fields and reflects the color-flavor locking property. In terms of the field E,-, the kinetic term of 
the Lagrangian (Q) then reads 
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/:kin = -^Tr 



We proceed by introducing both the vector field, = V^X°'/2, and the in-medium photon 
field, A^, by means of the covariant derivative 



D^^ = d,^-tgV,^-te^A,Q , (8) 

with Q = diag( 2/3, — 1/3, — 1/3 ) being the quark-charge matrix. The vector-meson field is 
rendered dynamical in the usual manner by introducing a kinetic term — (1/4) Tr F^^.F'^^, 
where F^i, = d^Vy — d^V^ + ig [V^, V^] is the field strength tensor. 

As noted in Ref. 0, the electromagnetic coupling in the CFL phase is modified. Since a 
diquark condensate carries baryon number, the physical (massless) photon field in medium is 
a linear combination of the original photon and the eighth component of the gluon field. The 
electric charge e in medium is therefore related to the vacuum charge e and strong coupling 



gs via e = e cos6 where cos^^ = y/S gs/ ^Sg^ + 4e^. Since, over a wide range of momenta, the 
strong coupling is much larger than the electromagnetic one, we will, for practical purposes, 
set cos ^ ~ 1 and drop the tilde notation for the photon field and the electromagnetic coupling 
in the following. The introduction of the electroweak interaction in the effective Lagrangian of 
the CFL and 2SC phases is presented in detail in Ref. |23j. 

The question now arises how to determine the the masses and coupling strengths of the 
vector mesons in the CFL environment. In Ref. ^T], where the HLS approach has first been 
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adopted to evaluate dilepton emission rates from CFL matter, vector-meson dominance (VMD) 
was assumed in connection with g = gp^^T, at its vacuum value (~ 6). From the KSRF rela- 
tion P3] it was then inferred that the vector-meson mass scales with the chemical potential ^q. 
On the other hand, in Ref. [12] a detailed weak-coupling calculation based on a Bethe-Salpeter 
equation has identified the axial/vector excitations (mesons) as diquark-hole states resulting 
in 

my = 2 A ^1 - e-^/s , (9) 

which for all practical purposes equals twice the gap (C = (3 — -\/3) vr^ ~ 30.7). In the 
recent work of Ref. [221 , (parametric) dependence of Eq.Q, my oc A, was corroborated 
by a nonperturbative argument, i.e., by matching higher-derivative contact terms of the low- 
energy effective theory, Eq. ((T)), with local terms obtained as a result of integrating out the 
heavy V^-fields in a generic chiral Lagrangian. Thus, we henceforth adopt my = 2 A also for 
the present case of a strong coupling environment, and do not distinguish between temporal 
and spatial vector- meson masses itlt and ms as was done in Ref. llj. The vector-meson mass 
term in the Lagrangian then acquires the standard vacuum form (my/2) V^V^. 

Within the matching procedure of Ref. |2S| it has furthermore been found that the vector- 
coupling strength scales according to^ 

This result suggests that, for large chemical potential, the self-interactions of the vector mesons 
are reduced with respect to (w.r.t.) the vacuum. By expanding the effective Lagrangian in 
the Goldstone boson fields, we see that g is identical to the coupling strength (7^7^77 of the 
vector meson to two pions. In particular, gp.,,T, scales in the same way as g. While in vacuum 
dprnr = g — G, the corresponding values in the CFL phase are parametrically suppressed. To 
obtain a more quantitative estimate for the pirn coupling, we follow Refs. [TTl in assuming 
that the KSRF relations remain valid in CFL matter, i.e., 

ml = 2ffl. (11) 

Here we have also assumed that VMD remains intact, implying that the HLS parameter has 
been fixed at a = 2.^ For the vector-meson interaction strength and their coupling to the 
Goldstone bosons we then have''' 

gp^^ = g^y^^. (12) 

^Using duality arguments, the analysis in '25' also exhibits that nucleon-type excitations are parametrically 
heavy with their masses proportional to /^/A. However, at moderate densities, where is not larger than the 
gap, this assertion no longer holds, and the (solitonic) "nucleon" or "delta" masses could become comparable 
to the vector meson ones; for large pairing gaps, it would thus be interesting to evaluate the effects of 7r-iV 
interactions - known to be strong in ordinary hadronic matter - in a CFL environment at temperatures above 
10 MeV. 

^In the weak-coupling analysis of Ref. fl^ it was found that the KSRF relations are modified at asymptotic 
densities, resulting in a=l. A similar feature arises within a renormalization-group treatment of the HLS 
approach at the finite temperature chiral phase transition |27| . 

^We would like to thank F. Sannino for pointing this out to us. 
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With a typical parameter choice of Aq = 50 - 150 MeV at /ig ~ 350 MeV (the latter implying 
/tt ~ 73 MeV), one obtains ^fpTTTr — 1 - 2.8, substantially smaller than in free space. 

For the reasons discussed above, we employ the vacuum form of the vector-meson mass 
term, implying that the vector-meson propagator is identical to the vacuum propagator. In 
order to preserve gauge invariance, it then turns out that also VMD needs to be applied in its 
vacuum version, i.e., without the distinction of temporal and spatial components in the 
vertex. Technically, this amounts to introducing into the Lagrangian a counter term of the 
form 

Aegf^ Tr QVo- vlA, Q Vi) - ^ e Tr (Q , (13) 

where the last term represents the standard V- 7 vertex^. Thus, after substituting the covariant 
derivative (jH)) in expression ((Tj) and expanding in the pion field vr, the complete interaction part 
of our Lagrangian for Goldstone bosons now reads 

CcBMt = -z^TT{Vo[n,do7r]+vm[7r,d.n]) + j{Al-vlA^)TT{[Q,7r]f 

Tr (Ao [Vo, vr] [Q, n] - vl A, [V,, tt] [Q, tt]) - ^ e A, Tr (Q V^) . (14) 

VMD reflects itself by the absence of a direct 7r7r7 interaction term. Finally, the physical 
photon and vector-meson field can be obtained by the usual diagonalization procedure. Note 
that, within the present approach, all parameters of the low-energy effective Lagrangian are 
now fixed. 



3 Photon Emission and Absorption 

Let us turn to the evaluation of the photon emission rates and mean free paths in CFL matter 
as a function of temperature and the value of the zero-temperature gap. As mentioned in the 
introduction, we restrict ourselves to processes within the SU{2) flavor sector, i.e., charged 
and neutral pions, as well as rho mesons. Since at temperatures of (9(1 MeV), e^e~ pairs 
have comparatively large number densities [TD^ , they will also be included in our analysis. For 
definiteness, we assume the validity of the BCS relation for the gap, A = Aq ^Jl — {T/T^"^, 
with a critical temperature Tc ~ 0.56 Aq. In order to demonstrate the dependence of our results 
on the (poorly known) gap, we perform our calculation with two different values, Aq = 50 MeV 
and Aq = 150 MeV. Unless otherwise stated, we set iiq = 350 MeV. 

3.1 Emission Rates and Emissivities 

The general set-up for the calculation of photon emission rates follows the ones performed 
in an ordinary hot meson gas in Refs. |2H1 1211 ■ The obvious differences expected in a CFL 
phase are due to the different masses and strong couplings of the pions and rho mesons, as 

similar procedure of restoring gauge invariance is often adopted when, e.g., including vertex form factors 
in hadronic effective theories '28!. The underlying microscopic origin of the counter terms remains open. 
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well as the modified pion dispersion relation. Both of these features lead to an increase in 
available phase space. The CFL pions are considerably lighter than in vacuum. Fixing the 
current quark masses at m„ = 3 MeV, = 7 MeV and = 120 MeV, we find m.^ ~ 12 
(38) MeV for Aq = 50 (150) MeV and temperatures of O (10 MeV). Throughout our numerical 
calculations, the temperature dependence of the pion mass, entering through A(T) in Eq. (jH), 
will be accounted for. 

The relevant processes which lead to one or two photons in the final state, are scattering of 
the form Mi M2 — > M3 7 (where Mj are pseudoscalar or vector mesons), annihilation. Mi M2 —>■ 
77, and rho-meson decays, p — > -KTcy. The (anomalous) neutral pion decay, 7r° ^77, has 
been considered in Ref. [9J for low temperatures of (9 (0.1 MeV), where the emissivity was found 
to be exponentially suppressed. For completeness, we have repeated the calculation for the 
larger temperatures we are interested in here. Note that purely electromagnetic processes, e.g., 
TT+TT^ — » 77, are relatively suppressed by a factor e'^/gp^^ ~ 10 - 100 as compared to strong 
processes such as tt+tt^ — > p° 7. At low temperatures, however, strong processes, requiring at 
least one in- or outgoing rho meson, are penalized by thermal suppression due to the relatively 
large rho mass. The in-medium Feynman rules and explicit expressions for the amplitudes 
are summarized in Appendix lA.ll and IA.21 respectively. We have verified that all amplitudes 
satisfy electromagnetic gauge invariance. 

In addition to the above processes, which can all occur in a vacuum environment, a novel 
photon source arises due to the medium-modified pion dispersion relation. With the in-medium 
pion velocity < 1, annihilation of the type 7r"'"7r~ — >■ 7 can become kinematically allowed, i.e., 
with all particles being on-mass shell. This is rendered possible due to the pion invariant-mass, 
p"^ = — (1 — f^)?^, being reduced and eventually turning spacelike at large 3-momentum. 
Thus, sufficiently fast pions can annihilate on slow ones. A more detailed discussion of this 
effect, which also proceeds in the inverse direction, thereby contributing to photon absorption, 
will be given in Sect. 13.31 

The photon emission rate for a given process, e.g., Mi{pi) M2{p2) — ^3(^3)7(9)5 is ob- 
tained from a phase space integral over the pertinent unpolarised scattering amplitude, A4, 
according to 

= W^^J J^WJ J2^^2J Wrm^^"^ Snipi+p2-ps-q) 



X 



M 



f{Ei,T)f{E2,T)[l + fiEs,T)], (15) 



where A/" is a (spin-isospin) degeneracy factor and f{Ei, T) are Bose-Einstein distribution 
functions. To evaluate the integral, we rewrite 

d^Pi f d Pi ^^^2 2 22 



(27r)3 2Ei J {2tt 



5{Et-mi-vtpi)QiEi), (16) 



with Vi being the velocity of particle i (different from 1 only for Goldstone bosons in the 
present work). One can then perform one of the 4-dimensional phase space integrations (for 
definiteness, i = 3) exploiting the (5-function in Eq. ()15|) . Since rotational symmetry persists 
in the medium, we may fix the photon momentum q in z-direction for subsequent purposes. 
The remaining angular dependences consist of the cosine of the polar angles, cos^i and cos 6*2, 
of particles 1 and 2 {w.r.t. the z-axis), and their relative azimuth = 0i — 02 (the integration 
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over the total azimuth $ = 0i + 02 becomes trivial) . With d^pj = \pj\'^ d\pj\d Qj and Ej dEj 
Vj \pj\ d \pj\ for j=l, 2, we arrive at 

dR 2-kM f \pi\dEidcoiiOi f \V2\dE2dcose2 ^, ^ , „ 
E— = / 2 / 2 "'^ + E2- E) 



d^q 8(27r)8y J vl 



X 



M 



2 

P3=pi+P2-q 

X /(El, T) /(E2, T) [1 + f{E, + E2-E, T)] . (17) 



6{{E, + E2- Ef - ml - vl (pi + P2 - q) 



For external (on-shell) particles, the energy integrations are restricted from below by Ei^^i-a = 
rrii. The integration can be performed using the remaining 5-function at its zero argument, 
(pQ, generating a Jacobian factor [2vf |pi| |p2| | sin^i sin 6^2 sin0o|]^^- The last four integrations 
are carried out numerically. 

The results for the emission rates in CFL matter for different gap values and temperatures 
are summarized in Fig. As a check of our numerical calculation, we have ensured agreement 
with the results of Refs. [28, 22| when using the vacuum values for masses (and velocities) of the 
pseudoscalar and vector mesons (pertinent curves for T = 50 MeV are displayed in the lower 
right panel of Fig. Q). Comparing the lower two panels, one observes that photon emission 
from CFL matter is strongly enhanced, which is essentially due to the increased phase space, 
as mentioned above. The particularly large rates for the processes 7r"'"7r~ — >■ 7 7 and np — >■ tt'j 
have their origin in the modified pion dispersion relation to be discussed in Sec. IH.HI The 
upper two panels indicate that strong processes involving rho mesons are further enhanced for 
smaller gap values, due to the reduced rho-meson mass. The emission rate for tt*^ ^ 7 7 is only 
shown in the vacuum case. In CFL matter this rate is kinematically suppressed: it vanishes 
if the photon energy is larger than 711^^/ ^Jl — v^. We finally remark that finite size effects of 
the effective meson states, i.e., vertex form factors, have not been included here (for a recent 
assessment of their suppression effects in photon production at high energies in normal hot 
hadronic matter, cf. Ref. [T^). 

The total photon emissivities, i.e., the emitted photon energies per unit volume and time, 
are now easily obtained by integrating the emission rate (jl5p (multiplied by the corresponding 
Bose enhancement factor) over the photon 3-momentum, 



Note that due to rotational symmetry the angular integration in the variable q can be done 
trivially yielding a factor 4 vr. The photon emissivity for all processes in CFL matter presently 
investigated is shown in Fig. |21 The process e+e~ ^77 dominates up to temperatures of 
about 3 (7) MeV for Aq = 50 (150) MeV. For larger T, 7i~^7[~ ^77 prevails. Strong 
processes involving rho mesons start to compete at temperatures of about 20 (25) MeV for 
Aq = 50 (150) MeV. However, the rho contributions are rather sensitive to the value of the rho- 
mass, subject to appreciable uncertainties due to, e.g., (i) sizable corrections to the leading- log 
estimate, Eq. Q, (ii) finite width effects (not included at present) from the p ^ nn coupling, 
which extend the p spectral strength down to the two-pion threshold. Therefore, our current 
calculations cannot yet reliably determine at which temperatures strong processes outshine the 
electromagnetic ones. Finally, let us point out that, due to the (partial) inversion of the mass 
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A,=^ 50 MeV, T=50 MeV Vacuum, T=50 MeV 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

E [GeV] E [GeV] 



Figure 1: Photon emission rates in CFL matter from various electromagnetic (thick hues) and 
strong-interaction (thin hues) sources as a function of the photon energy E. Upper plots: 
Ao = 50 MeV (left) and Aq = 150 MeV (right) at T = 20 MeV; lower plots: Aq = 150 MeV 
(left) and vacuum phase (right) at T = 50 MeV. 



hierarchy in the strange and nonstrange sector referred to in Sect. 2.1, we expect contributions 
from kaon pairs to be important or even dominant. These are not included at present, but can 
be treated by a straightforward extension of our calculations to SU{3)f- 



9 



Ao=50 MeV 



A,=1 50 MeV 




e e ->yy 
71 7t -> p Y 

71 p -> 71 Y 
p -> 71 71 Y 



10 15 
T [MeV] 



20 



10^ 



10' 



10^ 



o 

CO 



10' 



25 



10' 



10" 



'I 

''!■ 
'/' 

\l 



7/ 



11 



7t 7t -> Y Y 

7t"^7t"->Y 



n -> YY 



10 



15 
T [MeV] 



--- e e ->'Y'Y 
7t 7t -> p Y 

- - 7t p -> 71 Y 

- - - p -> 7t 7t Y 

20 25 30 



Figure 2: Photon emissivities in CFL matter from various electromagnetic (thick hues) and 
strong-interaction (thin hues) sources as a function of the temperature for Aq = 50 MeV (left) 
and Ao = 150 MeV (right). 



3.2 Mean free path 

To estimate the emitted photon luminosity of a finite size object an assessment of the pertinent 
mean free path is mandatory. In CFL matter, the latter is governed by Compton scattering 
off the light charged excitations, as well as photon annihilation reactions. In the SU{2) flavor 
sector, the lightest charged excitations are the pions as well as electrons and positrons. Similar 
to the emissivity discussed in the previous Section, at low T ~ 1 MeV one can expect pro- 
cesses involving to be dominant (due to their relatively large number densities), whereas 
towards higher T reactions involving pions ought to take over. As mentioned in the preceding 
subsection, another contribution arises from the novel "decay" process 7 tt^tt^, which owes 
its existence to the in-medium modified pion dispersion relation. 

The inverse mfp of a photon of energy E is given by its total cross section on all medium 
particles per unit volume (schematically (a^x^x) for each particle species X of density rix). 
For Compton scattering, •y^q) X'^^pi) — > 'y{q') X'^{p2), one has 

1 1 f d^pi f rf3p2 f d^q' 



ME) 2E I (2^0^ / (2^0^ / (2^ (2vr)^ S^^^ (g + Pi - g' - P2) 



X 



A^^^"-^^^" f{E,,T)[l±f{E2,T)][l + f{E',T)], (19) 



where the upper (lower) sign holds for X=tt (e) in connection with Bose-Einstein (Fermi- 
Dirac) distributions f{Ei,T) and f{E2,T). Here, the process amplitude is averaged over 
the initial photon polarisation and summed over all others. An analogous expression holds for 
77 e^e~ annihilation. 

In Fig. IHl we display our results for the mfp at a typical thermal energy of the incoming 
photon, E = 3T. We confirm that for Aq = 50 (150) MeV Compton scattering off prevails 
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Ao=50MeV Ao=150MeV 




5 10 15 20 25 5 10 15 20 25 30 

T [MeV] T [MeV] 



Figure 3: Photon mean free path from Compton scattering off tt^, e^, and photon annihilation 
into and tt^ pairs, for Aq = 50 MeV (left) and Aq = 150 MeV (right). 



for temperatures of up to about 5 (15) MeV, while at larger temperatures, charged pions are 
the more relevant scatterers. Photon "decay" into two pions can only occur if the photon 
energy E > 2m^/(l — v^), which for E = 3T implies a minimum temperature T^am = Tn.,^. 
Hence, this novel process does not play a significant role at low temperatures. As we will 
elaborate in the following Section, it becomes however relevant for T > m^r. This induces the 
rather sharp decrease of the solid line in Fig. IHl around T ~ 10 MeV (left) and T ~ 30 MeV 
(right), corresponding to the pion mass for Aq = 50 MeV and Aq = 150 MeV, respectively. 

3.3 In-Medium Pion Dispersion Effects 

We now return to the specific effects of the in-medium pion velocity which are at the origin of 
the novel emission and absorption processes, tt+tt^ — > 7 and 7 tt+tt^, respectively, already 
mentioned above. To illustrate their appearance in our calculations, let us for definiteness 
consider the reaction vr+p'' tt^'-j. This process has a contribution from a t-channel pion- 
exchange diagram. According to Eq. (|^. the pertinent propagator, D^, is modified by a 
factor v'l in the denominator, so that it takes the explicit form 

D-\Po, P) = Po - '"I p'-ml = {l-vl)-2E [E, - ^ El - ml cos Gi) . (20) 

As before, E is the final-state photon energy, Ei the initial-state pion energy, Bi the polar 
angle between their three momenta (the z-axis is chosen along the direction of the photon 
momentum) , and t = = {pi — qY is the four-momentum of the exchanged positively charged 
pion. In vacuum (f^ = 1), is always negative, i.e., in the spacelike regime. In the 
medium, however, the presence of the term E'^ (1 — v^) implies that can become timelike 
and therefore develop singularities in the physical region. Evidently, timelike values can only 
occur for sufficiently large photon energies. If, e.g.., the incoming pion is at rest, = 
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jt''" > — - •^/V^ Y y v/N^\;» - — >- 71 

p° — > — - - -> Tc"*" Ti"*" >--- -Krs^\y y 

Figure 4: Z-graph contribution to 7r+p° 7[~^^ and 7 vr"*" ^ 7 vr"*" . In time-ordered perturbation 
theory these diagram have a cut with three pion states. Arrows are pointing in the direction 
of momentum flow. 



requires a minimal photon energy of -Emin = 2m^/(l — v"^). Furthermore, in this situation 
the intermediate pion carries negative energy, readily interpreted as being its antiparticle with 
positive energy and reversed three momentum, that is, an incoming, on-mass-shell, vr". In the 
language of time-ordered perturbation theory this process represents a Z-graph contribution, 
depicted in the left panel of Fig. HI where the incoming 7r+ annihilates with an intermediate 
TT^ to produce the outgoing photon. The final scattering off the incoming p° to bring the vr"*" 
on-shell is no longer required, contrary to the vacuum case. In vacuum, the virtual vr"'" only 
goes on-mass shell in the limit of zero photon energy, in which case the phase space vanishes, 
and the emissivity receives no contribution. Completely analogous arguments apply for other 
charge combinations of the process tt p tcj. 

One can verify the above features more explicitly by directly analyzing the process vr(pi) + 
^(^2) lil)- Exploiting both energy and momentum conservation, Ei-\-E2 = E and P1-I-P2 = 
q, together with the on-mass shell dispersion relations for all three particles involved, one finds 
for cos 9i = ibl: 

E± = 3Ei±|pi|, (21) 

where we have used that = 1/3. This, in particular, implies that for an outgoing photon its 
(positive) energy E is always larger than Ei and therefore E2 > 0, that is, the Bremsstrahlungs- 
process n^pi) — > n{p2) 7(g) is not possible.^ 

A similar discussion of the in-medium pion dispersion effects pertains to processes that 
contribute to the mean free path. As we have already pointed out in the previous Section, 
at a temperature corresponding to the pion mass the pion-induced mean free path exhibits a 
rather sharp decrease. In our calculation this is associated with the w-channel diagram in 7 vr 
Compton scattering, which exhibits a behaviour analogous to the t-channel diagram discussed 
above, i.e., the intermediate pion becomes on-shell with negative energy. The interpretation 
via the corresponding Z-graph is that the incoming photon spontaneously decays into two 
on-shell charged pions, cf. the right panel of Fig. 0] 

For the numerical evaluation of the diagrams involving on-shell intermediate pions we in- 
clude in the (singular) pion-exchange propagators a finite width that can be attributed, e.g., to 
interactions with surrounding medium pions. A schematic estimate of the pion width is given 
in Appendix IA.3I For the present purpose it is sufficient to approximate the width by a con- 
stant; in the kinematical range of interest the evaluation of pTrvr interactions yields as an upper 
estimate a value of ~ 0.1 MeV for Aq = 50 MeV and F^ ~ 0.01 MeV for Aq = 150 MeV. 

^Equivalently, one can start from energy and momentum conservation for the Bremsstrahlungs-diagram to 
find that E2 must be negative, which again means that only the annihilation process can occur. 
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3.4 Comparison with Neutrino Properties 



At this point, a comparison with the results of Ref. JUl is in order, where an analogous 
investigation of neutrino emissivities and mean free paths was presented. On the one hand, 
we find that the photon emissivity exceeds the neutrino emissivity by more than ten orders 
of magnitude over the entire temperature region. On the other hand, as to be expected for 
weakly interacting neutrinos, their mean free path is larger than that of photons by several 
orders of magnitude. 

To confront the efficiencies of photon and neutrino emission in the radiation of energy from 
a hypothetical CFL star, we have evaluated the corresponding fluxes as follows. In order to gain 
a first estimate we ignore the possibility of a nuclear crust which might shield the CFL phase. 
For large opacities, as implied by the results found above, the photon fiux can be obtained 
from a convolution of the emissivity and mean free path, integrated over the photon energies. 
The neutrino fiux, for simplicity, has been estimated by multiplying the emissivities and mean 
free paths as given in the first entry of Ref. ^0], which amounts to the assumption that the 
neutrinos are thermalized when leaving the star.^° To facilitate a direct comparison with the 
findings of Ref. ^U], we have employed their values for the gap and the quark chemical potential, 
Ao = 100 MeV and fig = 400 MeV, respectively. The results are displayed in Fig. indicating 
that, despite of the extremely small mean free path of the photons, their emitted energy fiux 
exceeds that of neutrinos by one to three orders of magnitude over a wide temperature range. 

Note that the photon fiux approximately corresponds to that of a black-body emitter. This 
complies with the general expectation that, for very small photon mean free path, emission 
and absorption processes are close to equilibrium and refiect the surface temperature of the 
star. 



4 Early Thermal Evolution 

In this Section, we present a rough estimate of the early thermal evolution of a hypothetical 
CFL star. The thermal properties of CFL matter are expected to be dominated by the massless 
Goldstone boson 0, which results from the breaking of the baryon number symmetry U{1)b- 
The corresponding specific heat has been given in Ref. jU], 

= = 7.8 X 10^^ T^^^v erg cm'^ K'^ . (22) 

The thermal conductivity in a CFL phase has been estimated in Ref. [30] , where the contribu- 
tion from the massless Goldstone boson was found to be 

K^ = 1.2x 10^7 Tl,^y Acn, erg cm-^ s'' K'^ , (23) 

with Acm being the mean free path of the Goldstone bosons in cm. The dominant process 
for the mfp is the decay of Goldstone bosons into pairs of quark quasiparticles. Assuming 

^''This is a good approximation for temperatures above a few MeV [TOj . 

^^Here we are implicitly assuming that photons and neutrinos are emitted at the same temperature. This 
does not hold, however, once a significant temperature gradient has established at the surface. In that case, 
neutrinos will be predominantly emitted from hotter regions in the interior, due to their larger mean free path. 
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Figure 5: Comparison of photon and neutrino fluxes, with the parameter choice of Ref. [TUj . 
Ao = 100 MeV and ^lq = 400 MeV. 



thermalized Goldstone modes, the mean free path reads PU] 



A 



(T) 



4(21 - 8 In 2) 
15^2 ttT 



3A' 
2T 



(24) 



Note that A becomes large, i. e. greater than about 1 km, only for temperatures below 5 MeV. 

Since, at large T, the thermal conductivity is extremely small, we may estimate the cooling 
time of the CFL surface in the very early stage by neglecting the heat transport in the interior 
of the star. We thus use the thermal energy, t/th,sheii, stored in a shell with a thickness given 
by the photon mean free path, AR ~ A^, to write the luminosity as 



da 



th, shell 



dt 



(25) 



fAh.sheii is readily obtained by integrating the heat capacity over temperature and multiplying 
with the volume of the surface shell, 



rT 

Uth,she\\ = 4:71 R"^ A^ / c„,^(T') dT' . 
Jo 



(26) 



Using the black-body luminosity, = AtiR? aT^, together with Eq. (pUj) in expression (I2SI), 
and setting the radius of the star to i? = 10 km, we find that the cooling time At of the shell 
is of the order 



T- 

At ~ 10-^S In — , 



(27) 
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with Tj being the initial and Tj the final temperature. This cooling time scale is extremely 
short, and is to be compared with the diffusion time of the massless Goldstone bosons, cor- 
responding to the thickness AR of the shell. The latter time scale might be assessed by the 
simple relation (see, e.g., Ref. [12]) 



where tint — 10 s is the typical interaction time of the massless Goldstone bosons (estimated 
from their decay ^ gg) at T ~ 30MeV. With X^{T = 30MeV) ~ 10"^° cm we find 
t£) ~ 10~^'^s, which is much longer than the cooling time At. The large difference in these 
time scales justifies the neglect of the heat transport in the early stages of the cooling. 

5 Summary and Conclusions 

The main objective of this work has been an order of magnitude assessment of photon emission 
rates and mean free paths in a CFL phase of high-density quark matter. Our approach is 
based on a low-energy effective theory that describes the prevailing degrees of freedom at 
large densities and moderate temperatures, T ~ C (lOMeV). These are generalized Goldstone 
bosons that appear due to breaking of chiral and baryon number symmetries. We have followed 
the hidden local symmetry approach to introduce vector-meson excitations as well as the 
electromagnetic field in our analysis. The assumption of the KSRF relation in medium enabled 
us to relate the vector-meson self-interaction strength g and two-pion coupling gp^^^^ to the 
characteristic scales of the CFL phase, i.e., the gap A and the quark chemical potential fig, 
according to (^ptttt = g — 

V2A/f^, where U ^ 0.21 /i,. In the presence of in-medium modified 
pion dispersion relations, electromagnetic gauge invariance could be maintained by requiring 
both the p-7 vertex as well as the rho propagator (with mass rUp ~ 2 A) in its vacuum form 
(which emerged as a consequence of neglecting the distinction between temporal and spatial 
components of the rho mass). 

The resulting photon production and scattering amplitudes to lowest order have been convo- 
luted over available phase space in CFL matter using standard kinetic theory expressions. Our 
results suggest that photon emissivities might become very large, up to (9 (10^^ erg cm~^ s~^) 
at the highest temperatures expected to occur during (or immediately following) a supernova 
explosion, T ~ 30MeV [12^. Over a large range of CSC gap values, Aq = 50 - 150 MeV, we 
find that at temperatures below ~5 MeV the dominant process is e+e^ ^77, whereas for 
larger T electromagnetic pion annihilation, tt+tt" ^77, is most important. Processes invol- 
ving strongly interacting rho mesons could also play a major role. The temperature where the 
strong processes begin to outshine electromagnetic ones could not be accurately determined 
(although it might be well below 30 MeV), mostly due to uncertainties in the vector- meson 
mass spectrum. An interesting feature of the CFL phase is the appearance of the novel pro- 
cesses vr+vr" ^ 7, which are triggered by the in-medium pion dispersion relation due to the 
breaking of Lorentz invariance. 

We furthermore have found that the photon mean free path is extremely small over a wide 
range of temperatures. In particular, it is far below the typical radius of a compact star as 
long as the temperature is above about a few tens of keV. With a gap value of 50 (150) MeV, 
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the processes providing the largest opacity for T up to 5 (15) MeV are Compton scattering off 
thermal electrons and positrons, 76^ — 76^, as well as photon annihilation, 77 — e^e~, in 
agreement with previous calculations in a low-temperature approximation |3Tj. At larger T, 
Compton scattering off pions and photon annihilation into pion pairs are dominant. 

The comparison with analogous calculations of neutrino properties in CFL matter [TU] 
revealed that the photon emissivity is much larger than the neutrino one, whereas the mean 
free path is substantially shorter. A schematic estimate of the photon flux from a hypothetical 
CFL star with uniform temperature profile indicates that the photon flux could even exceed 
the volume-driven neutrino flux (temperature gradients will favor, however, the latter). 

We have also presented an order of magnitude analysis of the cooling time of the outermost 
surface layers of a (hypothetical) CFL star, in the very early stages of its evolution. Since, 
at the highest temperatures considered here, the thermal conductivity of CFL matter is very 
small, one may neglect heat transport in the interior of the star. The thus found cooling time 
is of the order 10~^^ s. This is to be compared with the diffusion time of a massless Goldstone 
boson, to ^ 10~^^ s, corresponding to the distance given by the photon mean free path at large 
temperature. These widely different time scales suggests that heat transport does not play a 
significant role in the very early cooling era. 

Let us finally comment on a few aspects of our work that have not been adequately treated 
and should be addressed in future analyses: 

(i) In our calculation of the emission rates and mean free paths we have restricted ourselves 
to the S'?7(2)-flavor case, i.e., did not include contributions from strangeness-carrying mesons 
in the pseudoscalar and vector octet. The analysis of Ref. ^H] showed that, while the K~ is 
slightly heavier than the pion, the K'^ is considerably lighter. Since the sum of the kaon masses 
is less than the sum of the pion masses, one can expect that the emission rates from K'^K" 
annihilation into photons will exceed the emission rates resulting from tc'^tc^ annihilation. 
Likewise, the kaon induced mean free path is expected to be smaller than the one induced 
by pions. Following our arguments of Sect. 13.41 we can assume, however, that the overall 
result for the flux will remain unchanged, being limited by black-body emission. At moderate 
densities, where the pion decay constant is of the order of the gap, one can also expect solitonic 
excitations in the baryon sector to become relevant. 

(ii) Our discussion of the cooling behavior of a compact star featuring CFL matter was 
limited to the very early stages of the thermal evolution, and we have furthermore made the 
simplifying assumption that the star accommodates a CFL phase which extends all the way 
up to the surface. We will postpone the investigation of the long-term cooling, including heat 
transport and more realistic scenarios in which the CFL phase is covered with a nuclear crust, 
to a future publication 122] • 
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A Appendix 

A.l In medium Feynman rules 

The Feynman rules follow from the effective Lagrangian discussed in Sect. |2] in a straight- 
forward fashion. From the kinetic term in Eq. (P) we obtain the pion propagator, 

^D. = - Y—-- (29) 

The vertex factors derive from the interaction part (fT^ . For the tt'^{p)p^ '^^{p') vertex 
we find 

ig/mTriP + p'),! , (30) 
where p = {po,vlp). The 7r+7r^p°7 contact vertex reads 

(31) 

and for 7r^7r'^p^7 we have 

(32) 

with g^^ = diag(l, —v^, —v^, —vl). The 7r+7r^77 contact term is 

-2ie^g^,y . (33) 

For the /0° - 7 coupling we obtain 

e m?/ , , 

-i—:^g,u. 34 

9 

Finally, the triple-rho vertex is given by 

i gpTTTT [gpv {p - p)x + gv\ {p - q)^t + gxp (g - p)^] , (35) 

where all momenta are flowing into the vertex. 
A. 2 Process amplitudes 

The process amplitudes for 7r+(pi) n^{p2) li.12) read 



^ ' [El- qiY - < (pi - qi)^ - 

M^,-^ ,. '^^'^r^rr-r (37) 

,pu 
'(3) 



Alfg" = 2 6^^'^". (3^ 
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By using crossing relations we immediately obtain the amplitudes for Compton scattering, 
7 7r^ — i> 7 7r^, and 77 tt+tt^, relevant for the calculation of the photon mean free path. 
The amplitudes for vr+(pi) vr^(p2) — ^ p'^iPs) lil) are 



(2pi-g)^(2p2-P3) 
{El - qf - vl (pi - q)2 - ml 

(2p2-g)M2pi-P3) 
{E2 - qY - vl (p2 - q)2 - ml 

A^JT) = '^egp^^r\ (41) 



■^(^ - "^''-^ - qy - vl (pi - q)2 - ' ^''^ 

_ {2p,-q)^i2pi-p,r .... 

-M(2) - e Qp^^ _ _ _ _ , i4Uj 



where the indices fi and refer to the photon and rho meson polarisation, respectively. The 
amplitudes for vr^p" 71^7 and p° 7r+7r^7 directly follow from the above expressions 
through crossing. 

For the processes '/r^(pi) vr°(p2) P'^iPs) lil) we find 



" ^^^^ {El - g)2 - ^2 (p^ _ q2) _ ^2 ' 



_ (2p3 + qT {Pi - P2)'^ - {Pi - P2Y (29 + p-iY + ^7^" (Pi - P2) ■ {q - Ps) 



(43) 

■^S) = -e^?p..^^^. (44) 

The amplitudes for vr^p^ — 71*^7, Tr'^p^ — > 7r^7 and p"*^ n'^n'^'y can again be obtained by 
crossing. 

A. 3 Pion width in a dense medium 

The width of the pion, generated by its interaction with surrounding medium pions, is defined 
by ImS^ = — mTrFTT. The pion self-energy can be obtained from the standard expression (see, 
e.g., Ref. |2S1) 

where f"^ is the Bose-Einstein distribution and J^'^'^'^'^'^ the forward-scattering amplitude of 
pions that scatter via (s-channel) rho-meson exchange. Upon evaluating Ai in the center-of- 
mass (cms) frame, Eq. takes the form 

S^(/co,k) = -— % / / {po,T) / — . (46) 



{2nY J Po "/smin |k|s-m2 + iImS 

The integration limits for the integral over the squared cms energy s are given by Smin,max = 
2ml + 2kopo±2 |k| |p|, and 

2 3 

ImSp = mprp = , q=--m^. (47) 



For simplicity, we here neglect effects of the in-medium pion dispersion relation. 
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